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The usual path integral formulation for scalar particles at finite density involves a sign problem, 
making numerical simulation impractical. We present alternative methods free of this difficulty. 
We apply these approaches to theory in 2 + 1 dimensions, presenting preliminary numerical 
results for physical quantities in the vicinity of the phase transition in the T — jj. plane. 
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1. Introduction 

A large class of fermionic and bosonic lattice theories exhibit severe sign problems at finite 
density, rendering numerical attempts to study such theories impractical. Attention has focused 
largely on fermionic theories and lattice quantum chromodynamics in particular, spurring numerous 
studies at chemical potentials /i < T. In this regime, the sign problem becomes less severe and can 
be dealt with by a variety of techniques. ^ At present, a generic solution to the sign problem that 
is valid in all regimes remains to be found. We take a less ambitious step toward this end, focusing 
on relativistic scalar field theories that exhibit sign problems. We sketch two related methods 
for solving the sign problem in the case of SO{N) models at finite temperature and density, and 
use these methods to study the phase diagram when N = 2 and in 2 + 1 dimensions. A detailed 
description of the methods, algorithms and results can be found in a future publication. 

We begin by considering a ?7(1) ~ S0{2) scalar field theory at finite chemical potential. Gen- 
eralization of our results to the global symmetry group SO{N) will be briefly discussed in Sec. ^. 
We work on an x j8 lattice (in lattice units, where a = 1) and impose periodic boundary con- 
ditions in the space and imaginary-time directions. The lattice action for this model, defined on 
Euclidean space-time, is given by 

= I t f^l-^-nP - rne-^^'''^n+.. " </'„*+ev^^'^ + V(|0„|)j . (1.1) 
n L V - 

The chemical potential pL has been introduced as a constant and imaginary zero-component of a 
gauge field, following the prescription of [^. At finite chemical potential, the action specified 
by Eq. |TTT] is complex. As a result, standard Monte Carlo techniques are inapplicable because 
the quantity no longer has an interpretation as a probability weight for generating statistical 

ensembles. We demonstrate how to avoid this issue by finding alternative representations for both 
the partition function and expectation values of operators in question. 

2. Alternative representations for Z(/i) 

We present two representations for the partition function, derived from strong coupling ex- 
pansion techniques, that avoid the sign problem. The first representation is based on the character 
expansion-an approach used in the more familiar context of lattice gauge theory. ^ After express- 
ing the complex scalar field in terms of radial and angular modes, one expands the weight in 
terms of the characters of [/(I). Subsequently, it is possible to integrate out the angular degrees of 
freedom. The resulting partition function takes the schematic form 

Zi(Ai)= £ [[pdp]e-'^'^P''^+<^^Q^^\ (2.1) 
{/n.v}eZ 







where the sum over integer- valued link fields I represents the sum over characters (Fourier modes) 
and p is the radial degree of freedom. The link field / represents the conserved current associated 



See [jl|] and references therein. 
For a review of these methods, see [ 
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with the U{1) symmetry of this model. As such, the sum over link fields is constrained to the 
subspace of physical configurations wfiich satisfy the divergence-free constraint 

£(V + /n,-v)=0, (2.2) 

V 

where Zn.-v = — 'n-ev,v- In this representation, the chemical potential couples to the charge Q, 
which is a functional of the link field I and given by 

i3e[/]=£/„,o. (2.3) 

n 

The action S\ is a complicated functional of the conserved current and radial mode. Nonetheless, 
one can verify that the action is real and therefore the partition function is given by a sum over real 
and positive weights, as desired for the purpose of Monte Carlo simulations. 

A second viable representation for the partition function is obtained from a hopping parameter 
expansion. Here, one Taylor expands the exponential of the nearest neighbor part of the action 
and then integrates term by term. Each term in the expansion is characterized by two link field 
configurations m and n. These link field configurations simply count the number of instances the 
nearest neighbor interaction and its complex conjugate appear in the expansion. Consequently, the 
Unk fields take on non-negative integer values. The partition function takes the form 

Z2(At) = L , (2.4) 

{m„.v}eN{n„,v}eN 

where the action 5^2 is real and, once again, the chemical potential couples to the total charge 
Q. In this representation, the charge is a functional of the difference m — n and thus we identify 
this quantity as the conserved current associated with the U{V) symmetry of our theory. One can 



explicitly verify that the difference m — n satisfies the divergence-free constraint given by Eq. 2.2. 

The methods described above easily generalize to the case of SO{N) theories with scalar 
fields transforming in the fundamental representation. Recall that the chemical potential is an anti- 
symmetric, real-valued matrix living in the adjoint representation of SO{N). With an appropriate 
change of basis, it is possible to block-diagonalize the chemical potential, where each block is a 
real 2x2 anti-symmetric matrix. If N is even, the theory may be regarded as a theory ofN/2 com- 
plex scalar fields coupled toN/2 independent chemical potentials and the aforementioned methods 
apply individually to each of the charged scalar fields. If N is odd, an additional real and neutral 
scalar field will be present. 



3. Algorithms 

We employ two algorithms for the purpose of simulating the partition functions Zi and Z2. The 
first algorithm is applied to Zi and is based on the Metropolis ^ accept/reject method. The main 
challenge of this approach is to perform updating of link variables in such a way as to maintain the 
constraint given by Eq. The problem is easily circumvented if one updates simultaneously link 
variables that form closed loops (conserved current loops) as shown in Fig. [l]; loops may either be 
local-forming the boundary of a plaquette, or global-forming a loop around the space-time torus. 
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Figure 1: Examples of global (left) and local (center) current loops. Arbitrary configurations (right) are 
obtained by multiple updates of local loops. 



The former is formally equivalent to solving the constraint by a change of variables to integer 
valued plaquettes [^. The chemical potential couples to fluctuations in the global current loops 
and consequently these fluctuations occur with a probability enhanced by a factor e^^ . Though we 
have not addressed the question here, it is possible to calculate expectation values of all neutral 
observables with this algorithm. 

We use the "worm" algorithm to perform numerical simulations of the partition function 
Z2; physical link field configurations are generated by the motion of a dynamical source (or sink) 
throughout the space-time lattice. The source traverses and consequently updates bond states with 
a probability governed by detailed balance. In this algorithm, only a subset of the configurations 
that are generated are physical; these correspond to closed loop configurations which arise when 
the source and sink occupy the same site. At finite chemical potential, the probability for forward 
motion of the source in the imaginary-time direction is enhanced by the factor . In a sequence 
of updates, it is possible and likely that the dynamical source will wind about the space-time torus 
resulting in an increase (or decrease) in the charge of the system at finite chemical potential. The 
worm algorithm provides a simple method for obtaining expectation values of link field dependent 
operators as well as the unnormalized distribution of 2-point correlation functions. While limited in 
utility in some respects, the algorithm is efficient in that it takes advantage of unphysical "fractional 
charge" intermediate states to undergo charge fluctuations, thereby avoiding tunneling processes 
that are inherent to the loop algorithm previously described. 



4. An application 

We apply the above formalism to a 2 + 1 dimensional complex scalar field theory at finite 
temperature and density with the potential V = m^\^\^ + Numerical simulations of this 

model are performed using both the loop and worm algorithms with the aim of understanding the 
phase diagram in the T — jj. plane and in the vicinity of the phase transition. We present preliminary 
results of the effort. 
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Figure 2: f/ as a function of m^ for j3 ,L = 20 — 36. 
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Figure 3: n//X^ as a function of j3,Lfor/i =0.1 -0.5. 



At zero temperature and density, the theory under consideration is expected to exhibit a second 
order phase transition, located by first fixing the self-coupling at A = 192 and then tuning the mass 
parameter to the critical value m^. The critical value for the mass parameter is determined by 
studying the Binder cumulant 



1 



HQ- 



2\2 



(4.1) 



as a function of m^. At large finite volumes, the charge distribution is Gaussian when m < and 
exponential when m > nic. Consequently, in the thermodynamic limit U (m^) will approach zero in 
the broken phase and diverge in the unbroken phase. Results obtained via the worm algorithm are 
displayed in Fig. ^ for linear space-time volumes j8 , L = 20 — 36. We estimate the critical mass from 
the intersection of level curves at large space-time volumes, finding f» —26.05 for A = 192. 

Near the critical point (m^, A) we investigate the dependence of the number density as a func- 
tion of chemical potential. According to [^, the number density should depend quadratically on the 
chemical potential and the curvature ^ is expected to be a universal and nonperturbative constant. 
We determine the curvature by plotting the ratio n/}X^, where n is the number density, as a function 
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Figure 4: n//i as a function of temperature for = 0. 1 — 0.5. 




Figure 5: j as a function of jU for L = 12 — 36. 



of the linear space-time volume for values of the chemical potential ranging from pi =0.1 — 0.5. 
Results shown in Fig. || were once more obtained using the worm algorithm. For pL > 0. 1, the level 
curves converge to .'^ k, 0.32 in apparent agreement with an indirect determination of ^ = 0.32(1) 
by [i- 

Using the loop algorithm, we study the qualitative behavior of the number density as a function 
of temperature for chemical potentials ranging from = 0. 1 — 0.5. At temperatures T <C /X, the ef- 
fective theory for the Goldstone mode obtained by the techniques of [^] predicts weak dependence 
of the number density on temperature. Specifically, one finds 



n{pL) = mpL^ i + e^{T/pLf 



(4.2) 



which appears in agreement with our numerical calculations displayed in Fig. ^ Furthermore, we 
expect the shift in the critical temperature to depend linearly on the chemical potential JTOj], where 
the constant of proportionality ^ is universal and nonperturbative. We can determine the value 
of ^ by studying the chemical potential dependence of the susceptibility associated with |0p for 
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fixed j8. The susceptibility is defined by 

;t = /3L^((WV(W)(W)) , (4.3) 

where |(^p is the space-time volume average of An estimate of is attainable from the 
location of the maxima of the susceptibihty /Xm(^) as one approaches the thermodynamic limit 
L —>■ oo. In this limit, the susceptibility is expected to diverge at the phase transition and therefore 
we expect "^^^ = j3/im(°°)- Preliminary results for the susceptibility at linear spatial volumes 
L = 12 — 36 and inverse temperature j8 = 6 are shown in Fig. |5[ A precise determination of ^ is in 
progress. 

5. Conclusion 

We have presented two alternative representations for the partition function for scalar field 
theories at finite density and temperature. These representations are free of sign problems and are 
applicable to the numerical study of SO{N) models with fundamental matter. We have outlined 
several algorithms for simulating these partition functions and have applied the methods to 
theory in 2 + 1 dimensions. Preliminary numerical results for the universal constant ^ appear in 
agreement with previous studies and our finite temperature results are qualitatively consistent with 
analytic calculation. 

It is unlikely that the expansion techniques we have discussed are applicable to fermionic 
theories exhibiting sign problems. Nonetheless, we believe that the pursuit of other representations 
for fermionic theories is worth consideration. 
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